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We have numerically calculated the single band Wannier functions for interacting Bose gases in 
optical lattices with a self-consistent approach. We find that the Wannier function is broadened by 
repulsive atom interaction. The tunneling parameter J and on-site interaction U computed with 
the broadened Wannier functions are found to change significantly for different atomic number per 
site. Our theory can explain the nonuniform atomic clock shift observed in [Campbell et al, Science 
313, 649 (2006)]. 
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I. INTRODUCTION 

In condensed matter physics, people are often inter¬ 
ested in the ground state of the system and its low en¬ 
ergy excitations. This allows people to focus only on the 
lowest bands of the systems by mapping the system to a 
lattice model with Wannier functions [T]. The Hubbard 
model is arguably the most famous of all such lattice 
models [2]. In addition, Wannier functions are a more 
natural and better choice than Bloch waves for narrow- 
band materials in computational physics [3]. There is a 
lot of freedom to choose a set of Wannier functions as ba¬ 
sis and there has been great efforts to find the best Wan¬ 
nier functions [4]. However, all the discussion is done 
in the context of single-particle physics. The broaden¬ 
ing of Wannier functions by repulsive interaction seems 
have never been discussed in traditional condensed mat¬ 
ter physics. 

The situation starts to change with the development 
of ultracold atomic physics, in particular, with the ex¬ 
perimental observation of superfluid to Mott transition 
in optical lattice EE]. In the early theoretical treat¬ 
ment of such a system, the mapping from the realistic 
continuous system to the Bose-Hubbard model is done 
with single particle Wannier function [ 7 ]. However, due 
to the simplicity of the periodic potential and small en¬ 
ergy scales in the system, it is tempting to think that 
the broadening of Wannier function by the on-site repul¬ 
sive interaction may have significant effect on the system. 
There have been many theoretical efforts [8IIT3] trying to 
give a good description of interaction effects on Wannier 
function. 

There are also strong experimental evidences on the 
the broadening of Wannier functions. In atomic clock 
experiment pTHlT] , atomic interaction is the main reason 
for frequency shift. In the experiment of atomic 

clock in optical lattice m, nonuniform frequency shift 
was observed for different occupation number per site. 
This is clearly due to the broadening of Wannier function 
by the repulsive atomic interaction. 

In Ref. [18] , a self-consistent approach is developed to 
take account of the interaction effect on Wannier func¬ 


tion. The use of a different set of Wannier functions will 
result in a different tunneling parameter J and the on-site 
interaction U for the lattice model and thus a different 
ground state. The self-consistent approach in Ref. m 
uses a general variational principle to choose the set of 
Wannier functions that minimizes the ground state en¬ 
ergy of the lattice model. This method is in spirit the 
same as the MCTDHB theory [TM2T] . 

In this work we use the self-consistent approach in 
Ref. HHI to compute the interaction broadened Wan¬ 
nier functions for a Bose gas in an optical lattice. We 
focus on both the superfiuid regime and the Mott insula¬ 
tor regime. The Wannier functions are used to calculate 
the tunneling parameter J and the on-site interaction U 
in the Bose-Hubbard model. They are found to be signif¬ 
icantly affected by the s-wave scattering length, lattice 
strength, and most importantly the number of particles 
per site. In the end, we apply the approach to the ex¬ 
periment in Ref. [T5| ; our theoretical results match very 
well the experimental data. 

Our paper is organized as follows. In section H, a quick 
review of the self-consistent approach is given in Ref. 
m- In section HI the nonlinear equations are solved 
in the superfluid regime for one dimensional optical lat¬ 
tice; the tunneling parameter J and on-site interaction 
U are calculated with the changing of lattice depth and 
interaction strength. In section IV, we solve the nonlin¬ 
ear equations for the Mott insulator regime where the 
expansion basis are single particle Wannier functions; J 
and U are calculated accordingly. In section V, the the¬ 
ory is applied to the experiment in Ref. m and a good 
agreement is found between theory and experiment. 


II. SELF-CONSISTENT APPROACH FOR 
WANNIER FUNCTIONS 

In this section we give a brief summary on the self- 
consistent approach to compute the interaction effects 
on Wannier functions developed in Ref. [18] . We consider 
a Bose gas where the weak atomic interaction can be well 
described by the 5-wave scattering. The second quantized 
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Hamiltonian for this kind of system is given by 


H = 


[ dr'ipUr)[—-^\/‘^-\-V(r)]'ip(r) 

J 2m 


' 2 


J dr[ip\r)'ip\r)'ip{r)'ip{r )], 


( 1 ) 


where m is the atomic mass, V{r) is the external poten¬ 
tial. We choose V (r) = Vosin{7rr)‘^ in this context, where 
Vo is the optical lattice strength. And go = as jm is 
the interaction strength related to the s-wave scattering 
length ttg. The single-band approximation is to expand 
the bosonic field operator 'ip{r) as 


^{r) = Y,ajWj{r), (2) 

j 


where Wj{r) = W{r — rj) is the first band Wannier func¬ 
tion at site j and dj is the associated annihilation opera¬ 
tor. The ground state \Gt) in the single-band approxima¬ 
tion can be generally written as \Gt) = F(dt)|vaccum), 
where F is the function to be found by solving the re¬ 
sulted lattice model. The ground state energy Eg = 
{Gt\H\Gt) certainly changes with the choice of Wannier 
function Wj{r). The best Wannier function is the one 
that minimize the single-band ground state energy Eq. 
Mathematically, this is to do the following variation 

6Eg _ . . 

5W^{r) 6W*{r) “ ’ ^ ^ 

with the orthonormal constrains 

hj = J drW*{r)W{r - rj) = 5oj • (4) 

/ij’s are the usual Lagrangian multipliers. With straight¬ 
forward computation, a nonlinear equation was obtained 
for the interacting Wannier functions m 


HW{r -rj)=Y, {a],aj,)HoW{r + Vj, - 

j jlj2 

JsjA 

+ 90 

jlj2 

xW{r + - VjOW (r + - Vj^). 



FIG. 1: Comparison of the ground state energies of Bose- 
Hubbard model between the usual single-particle method and 
our self-consistent method. The self-consistent method pro¬ 
duces a lower ground state energy. Er = h^k\/2m is the 
recoil energy, /cl is the wave vector of optical lattice laser. 


and 

U = go J dr\W{r)\\ (8) 

Eq. <§ and Eq.(|^ need to be solved self-consistently to¬ 
gether to find the best Wannier function W{r). 



FIG. 2: Blue dashed line is single particle Wannier function, 
and red line is interacting Wannier function by self-consistent 
method, x is in unit of A/2, where A is the wave length of 
optical lattice laser. 


( 5 ) 

where (•) represents averaging over the ground state of 
the system. The ground state is to be found with the 
Bose-Hubbard model 

-ff?. =+ ^y]a1'ai(a|ai - 1), (6) 

(b'> « 

where 

J = - [ drW^r - rj)HoW{r - r^_i) (7) 


As an example, we consider a one dimensional six-site 
Bose-Hubbard model with one atom per site so that we 
can find its ground state with exact diagnahzation[22] . 
We use two different ways to compute J and U in the 
model: (1) with the single particle Wannier function; 
(2) with the interactin g Wa nnier function obtained self- 
consistently with Eqs.( |5|6| ). The ground state energy 
of this Bose-Hubbard model is compared for these two 
methods in Fig. The energy computed with the self- 
consistent method is indeed lower. Fig. shows one 
Wannier function that we obtained in superfiuid regime. 
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which is apparently broadened due to the interaction. 
These broadened Wannier functions can influence the 
tunneling parameter J and on-site interaction U in the 
single band Bose-Hubbard model. In the following two 
sections, we shall compute the broadened Wannier func¬ 
tion in both superfluid regime and Mott regime for one 
dimensional systems. 

III. SUPERFLUID REGIME 

In the superfluid regime all the particles are condensed 
into the ground state of the system which is a Bloch 
state; it is more convenient to use the Bloch basis. The 
nonlinear equations for Bloch waves are m 

i'k'fpkir) = Hoipk{r) 

+ 90 Pk^kksk^'ip*k,{r)'ipk3{r)'ipkAr), ( 9 ) 

{kikksk^) 

where Pkikkskt = {Vk^k^kibk^)/{VJ)k) with bk = 
‘ In the superfluid phase, the Bogoli- 
ubov mean-field theory [23] can be used to determine and 
compute Pkikk 3 k 4 and other coefficients in Eq.(|^ [TS] . 
There is no indeed to solve Eq. 

In the computation, we expand the Bloch function 
with plane waves, 

Mr) = + (10) 

where U is the volume of a cell and N is the number of 
cells. After plugging the above equation into Eq.(|^, we 
will obtain a set of nonlinear equations for a{k). We solve 
these nonlinear equations numerically and then construct 
the Wannier functions of different energy bands by care¬ 
fully choosing the phases of all the Bloch functions with 
Kohn’s method [24]. 

Eig. shows that J and U change with lattice strength 
Vb and interaction parameter go for different mean par¬ 
ticle number on a single site no- It is apparent that the 
atomic number per site does not change the overall trends 
but does change J and U significantly. As the result 
of broadening of Wannier functions, J increases with no 
while U decreases. In the range that we show in Eig. 

J with our approach can change up to 100% compared 
to that of a single particle result at no = 3 while U can 
change up to 30%. 

IV. MOTT INSULATOR REGIME 

In deep Mott-insulator regime the ground state can 
be approximated with |no, no,..., no), where no is atom 
number per cite. As a result, the nonlinear equation 
for interacting Wannier function in the Mott regime is 
simplified to 


E = ^oW"(r) + 5o(no - l)\W{MW{r) 

j 

+ 2gono E mr-rj)fW{r). 

Vj^O 

( 11 ) 

We expand the Wannier function in terms of single¬ 
particle Wannier functions on the same site and its near¬ 
est neighbors, 

M 

W{r-rj)=Y,[cnMr-rj-i) 

n=l 

+ bnWn{r - Tj) -h CnWn{r - + 

where Wn{r — Vj) is the single particle Wannier function 
for band n and site j. In this numerical calculation, we 
set M = 3 and dimension for this section is ID. With 
this expansion, we solve Eq.([TT]) numerically to find the 
interaction-broadened Wannier functions, and compute J 
and U. The results are shown in Fig. It is clear that 
the general trends that J and U change with go and Vb in 
the Mott regime are similar to the ones in the superfluid 
regime. However, there are differences. Specifically, as 
shown in Fig. |^a,c), the change of both J and U with no 
has little dependence on the lattice depth Vb- As shown 
in Fig. J can change up to 32% in our self-consistent 
approach at no = 4 compared to that of single particle 
Wannier function while U can change up to 14%. 

V. EXPLANATION OF NON-UNIFORM CLOCK 
SHIFT 

In atomic clock measurement, the clock frequency can 
shift due to collision of atoms. In an experiment reported 
in Ref. [15], the atomic clock shift of Rh was measured. 
In the experiment, a Rh Bose-Einstein condensate was 
prepared in the |F = l,m/ = —1) state and loaded into a 
3D optical lattice. With the increase of the lattice depth, 
the system changed from superfluid phase to Mott insu¬ 
lator (MI) phase. Due to the trapping potential, the 
atomic gases was separated into MI shells each of which 
has a different occupation number no per site. One can 
use radio wave to excite atoms in F = 1 state to F = 2 
state. In different hyperflne states, the scattering lengths 
between atoms are different. Therefore, the atoms trans¬ 
ferred to the F = 2 state have a slightly different mean 
field energy; this can cause a clock frequency shift, which 
is given by [23l [25| 

5v =P^{a 2 i-aii)/aii, (13) 

where an and a 22 are scattering lengths for atoms in 
F = 1 and F = 2 state, respectively, and ai 2 is scat¬ 
tering length between atom in F = 1 state and atom in 
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FIG. 3: Superfluid regime, (a) Change of tunneling parameter J with potential depth for different no; (b) Change of J with 
interaction strength for different no; (c) Change of on-site interaction U with potential depth for different no; (d) Change of U 
with interaction strength for different no. 





FIG. 4: Mott insulator regime, (a) Change of tunneling parameter J with potential depth for different no; (b) Change of J 
with interaction strength for different no; (c) Change of on-site interaction U with potential depth for different no; (d) Change 
of U with interaction strength for different no. 
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FIG. 5: Clock shifts for different atomic occupation numbers 
per site. The shift becomes smaller for higher occupation. 

F = 2 state. If the on-site interaction U is calculated 
with single-particle Wannier functions, this clock shift is 
independent of no, the number of atoms per site. How¬ 
ever, it was observed in the experiment, the clock shift 
decreases with no as shown in Fig. 

In our view, this decrease of clock frequency shift in 
Fig. 1^ is caused by the broadening of Wannier func¬ 
tion: when no increases, the Wannier function becomes 
broader and U decreases as shown in Fig. [^c,d); con¬ 
sequently, the clock shift decreases. With our self- 
consistent method, we calculated the 3D interacting 


Wannier function for different atom number no per site 
in 3D Mott-insulator regime, and eventually U and the 
clock shift 6u. The results are compared to the experi¬ 
mental results in Fig. there is a very good agreement. 


VI. CONCLUSION 

With a self-consistent theory developed earlier, we 
have computed the effect of interaction on single-band 
Wannier function. In particular, we have considered both 
superfluid regime and Mott regime. We found that as 
the result of broadening of Wannier function by interac¬ 
tion, the tunneling parameter J and U can change sig¬ 
nificantly. Our theory was applied to a clock shift ex¬ 
periment: a very good agreement was found between our 
theoretical results and experimental results. The regime 
near the transition from superfluid to Mott insulator is 
not studied in this work and will be investigated in the 
future. 
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